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Abstract 

We develop the shape derivative analysis of solutions to the problem of scattering 
of time-harmonic electromagnetic waves by a bounded penetrable obstacle. Since 
boundary integral equations are a classical tool to solve electromagnetic scattering 
problems, we study the shape differentiability properties of the standard electro- 
magnetic boundary integral operators. To this end, we start with the Gateaux 
differentiability analysis with respect to deformations of the obstacle of boundary 
integral operators with pseudo-homogeneous kernels acting between Sobolev spaces. 
The boundary integral operators of electromagnetism are typically bounded on the 
space of tangential vector fields of mixed regularity TH 5 (divr, T). Using Helmholtz 
decomposition, we can base their analysis on the study of scalar integral operators 
in standard Sobolev spaces, but we then have to study the Gateaux differentiability 
of surface differential operators. We prove that the electromagnetic boundary in- 
tegral operators are infinitely differentiable without loss of regularity and that the 
solutions of the scattering problem are infinitely shape differentiable away from the 
boundary of the obstacle, whereas their derivatives lose regularity on the boundary. 
We also give a characterization of the first shape derivative as a solution of a new 
electromagnetic scattering problem. 

Keywords : Maxwell's equations, boundary integral operators, surface differential 
operators, shape derivatives, Helmholtz decomposition. 
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Introduction 



Consider the scattering of time-harmonic electromagnetic waves by a bounded obstacle 
Q in M 3 with a smooth and simply connected boundary T filled with an homogeneous 
dielectric material. This problem is described by the system of Maxwell's equations, 
valid in the sense of distributions, with two transmission conditions on the boundary or 
the obstacle guaranteeing the continuity of the tangential components of the electric and 
magnetic fields across the interface. The transmission problem is completed by the Silver- 
Miiller radiation condition at infinity (see |23j and [24]). Boundary integral equations are 
an efficient method to solve such problems for low and high frequencies. The dielectric 
scattering problem is usually reduced to a system of two boundary integral equations for 
two unknown tangential vector fields on the interface (see [6] and |24]). We refer to [9] and 
|10| for methods developed by the authors to solve this problem using a single boundary 
integral equation. 

Optimal shape design with the modulus of the far field pattern of the dielectric scat- 
tering problem as goal is of practical interest in some important fields of applied mathe- 
matics, as for example telecommunication systems and radars. The utilization of shape 
optimization methods requires the analysis of the dependency of the solution on the shape 
of the dielectric scatterer. An explicit form of the shape derivatives is required in view 
of their implementation in a shape optimization algorithms such as gradient methods or 
Newton's method. 

In this paper, we present a complete analysis of the shape differentiability of the 
solution of the dielectric scattering problem using an integral representation. Even if nu- 
merous works exist on the calculus of shape variations |14[ [25j [26j l3T| [32] , in the framework 
of boundary integral equations the scientific literature is not extensive. However, one can 
cite the papers [27], [29 J and [28], where R. Potthast has considered the question, starting 
with his PhD thesis |30| . for the Helmholtz equation with Dirichlet or Neumann boundary 
conditions and the perfect conductor problem, in spaces of continuous and Holder con- 
tinuous functions. Using the integral representation of the solution, one is lead to study 
the Gateaux differentiability of boundary integral operators and potential operators with 
weakly and strongly singular kernels. 

The natural space of distributions (energy space) which occurs in the electromagnetic 
potential theory is TH ^ (divr, T), the set of tangential vector fields whose components 
are in the Sobolev space H~z (T) and whose surface divergence is in H~ 2 (r). We face two 
main difficulties: On one hand, to be able to construct shape derivatives of the solution - 
which is given in terms of products of boundary integral operators and their inverses - it is 
imperative to prove that the derivatives are bounded operators between the same spaces 
as the boundary integral operators themselves. On the other hand, the very definition of 
shape differentiability of operators defined on TH~2 (divr, T) poses non-trivial problems. 
Our approach consists in using the Helmholtz decomposition of this Hilbert space. In this 
way, we split the analysis in two main steps: First the Gateaux differentiability analysis 
of scalar boundary integral operators and potential operators with pseudo-homogeneous 
kernels, and second the study of derivatives with respect to smooth deformations of the 
obstacle of surface differential operators in the classical Sobolev spaces. 
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This work contains results from the thesis |21j where this analysis has been used to 
develop a shape optimization algorithm of dielectric lenses in order to obtain a prescribed 
radiation pattern. 

The paper is organized as follows: 

In section Q] we define the scattering problem of time-harmonic electromagnetic waves 
at a dielectric interface and the appropriate spaces. In section [2] we recall some results 
about trace mappings and boundary integral operators in electromagnetism, following the 
notation of [U [23]. We then give an integral representation of the solution following |9j. 
In section (3[ we introduce the notion of shape derivative and its connection to Gateaux 
derivatives. We also recall elementary results about differentiability in Frechet spaces. 

The section^] is dedicated to the Gateaux differentiability analysis of a class of bound- 
ary integral operators with respect to deformations of the boundary. We generalize the 
results proved in [27, [29] for the standard acoustic boundary integral operators, to the 
class of integral operators with pseudo-homogenous kernels. We also give higher order 
Gateaux derivatives of coefficient functions such as the Jacobian of the change of variables 
associated with the deformation, or the components of the unit normal vector. These re- 
sults are new and allow us to obtain explicit forms of the derivatives of the integral 
operators. 

The last section contains the main results of this paper: the shape differentiability 
properties of the solution of the dielectric scattering problem. We begin by discussing the 
difficulties of defining the shape dependency of operators defined on the shape-dependent 
space TH 5 (divr, T), and we present an altervative using the Helmholtz decomposition 
(see [11] ) on the boundary of smooth domains. We then analyze the differentiability of 
a family of surface differential operators. Again we prove their infinite Gateaux differen- 
tiability and give an explicit expression of their derivatives. These results are new and 
important for the numerical implementation of the shape derivatives. Using the chain 
rule, we deduce the infinite shape differentiability of the solution of the scattering prob- 
lem away from the boundary and an expression of the shape derivatives. More precisely, 
we prove that the boundary integral operators are infinitely Gateaux differentiable with- 
out loss of regularity, whereas previous results allowed such a loss [28] , and we prove that 
the shape derivatives of the potentials are smooth far from the boundary but they lose 
regularity in the neighborhood of the boundary. 

These new results generalize existing results: In the acoustic case, using the variational 
formulation, a characterization of the first Gateaux derivative was given by A. Kirsch in 
[20] for the Dirichlet problem and then for a transmission problem by F. Hettlich in 
[15} 116] . R. Potthast used the integral equation method to obtain a characterisation of 
the first shape derivative of the solution of the perfect conductor scattering problem. 

We end the paper by formulating a characterization of the first shape derivative as 
the solution of a new electromagnetic scattering problem. We show that both by directly 
deriving the boundary values and by using the integral representation of the solution, we 
obtain the same characterization. 
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1 The dielectric scattering problem 



Let Q denote a bounded domain in R 3 and let Q c denote the exterior domain R 3 \f2. In 
this paper, we will assume that the boundary T of is a smooth and simply connected 
closed surface, so that f2 is diffeomorphic to a ball. Let n denote the outer unit normal 
vector on the boundary T. 

In f2 (resp. f2 c ) the electric permittivity e« (resp. e e ) and the magnetic permeability 
Hi (resp. fj, e ) are positive constants. The frequency u is the same in Q and in Q c . The 
interior wave number «j and the exterior wave number n e are complex constants of non 
negative imaginary part. 

Notation: For a domain G C t 3 we denote by H S (G) the usual L 2 -based Sobolev 
space of order s G R, and by Hf QC (G) the space of functions whose restrictions to any 
bounded subdomain B of G belong to H S {B). Spaces of vector functions will be denoted 
by boldface letters, thus 

W{G) = {H s {G)f . 
If D is a differential operator, we write: 

H S (D, n) = {u G H s (n) : Du G H s (ft)} 
Hf oc (D,rF) = {u G Hf oc (JF) : Du G Hf oc (JF)} 

The space H s (D,f2) is endowed with the natural graph norm. When s = 0, this defines 
in particular the Hilbert spaces H(curl, $7) and H (curl curl, Q). We denote the L 2 scalar 
product on T by (•, -}r- 

The time-harmonic Maxwell's sytem can be reduced to second order equations for the 
electric field only. The time-harmonic dielectric scattering problem is then formulated as 
follows. 



The dielectric scattering problem : Given an incident field E mc G Hi oc (curl,] 

that satisfies curl curl E mc — K 2 E mc = in a neighborhood of Q, we seek two fields 
E l G H(curl,Sl) and E s G Hi oc (curl, S7 C ) satisfying the time-harmonic Maxwell equations 



curl curl E* 
curl curl E s 








in 0,, 
in ft c , 



(1.1) 
(1.2) 



the two transmission conditions, 



nxE' = nx(E s | E inc ) on T 

H~\n x curl E*) = x curl(E s + E inc ) on T 



(1.3) 
(1.4) 



and the Silver-Muller radiation condition: 



lim \x\ 



curlE s (x) 



x 

X 



iK e E s (x) 



(1.5) 
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The interior and exterior magnetic fields are then given by H* = E* and H s = E s . 

It is well known that this problem admit a unique solution for any positive real values of 
the exterior wave number [6l [2TJ [23] . 

An important quantity, which is of interest in many shape optimization problems, is 
the far field pattern of the electric solution, defined on the unit sphere of M 3 , by 

E°°(x) = lim 47r|x| - ^ , with ^— = x. 

2 Boundary integral operators and main properties 
2.1 Traces and tangential differential calculus 

We use surface differential operators and traces. More details can be found in [8| 124] . 

For a vector function v £ (^ fc (M 3 )) <? with k,q £ N*, we note [Vv] the matrix whose 
the i-th column is the gradient of the i-th component of v and we set [Dv] = T [Vv]. The 
tangential gradient of any scalar function u £ ^(T) is defined by 

V r n = VCt| r - (Vu| r • n) n, (2.1) 

and the tangential vector curl by 

curlr u = Vn|r x n, (2.2) 

where u is an extension of u to the whole space M 3 . For a vector function u £ (^(r)) 3 , 
we note [Vru] the matrix whose the i-th column is the tangential gradient of the i-th 
component of u and we set [Dpu] = r [Vru]. 

We define the surface divergence of any vectorial function u £ ( < ^ ?fc (r)) 3 by 

divr u = div Q|p — ([Viim]n • n) , (2.3) 

and the surface scalar curl curlp r u r by 

curlp u = n • (curl u)) 

where u is an extension of u to the whole space M 3 . These definitions do not depend on 
the extension. 

Definition 2.1 For a vector function v £ ( c rf° c (Q,)) 3 and a scalar function v £ < ^°°(f2) 
we define the traces : 

iv = v\ r , 
7£>v := (n x v)| r (Dirichlet) and 
7at k v := K _1 (n x curlv)| r (Neumann). 
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We introduce the Hilbert spaces H S (T) = 7 ^H s+ ^(Q)^j , and TH S (T) = 7D (h s+ ^(0)^ 
For s > 0, the traces 

7 : H s+ ^{n) -> # s (r), 

7 D : H s+ 5(0) -> TH s (r) 

are then continuous. The dual of H S (F) and TH s (r) with respect to the L 2 (or L 2 ) scalar 
product is denoted by H~ S (T) and TH~ s (r), respectively. 

The surface differential operators defined here above can be extended to the Sobolev 
spaces: The tangential gradient and the tangential vector curl are linear and continuous 
from H S+1 (T) to TH s (r), the surface divergence and the surface scalar curl are linear 
and continuous from TH S+1 (T) to H S (T). 

Definition 2.2 We define the Hilbert space 

Tl-H(divr,r) = jj G TH-2(r),div r j G H~^(T)\ 

endowed with the norm 



11 m th~2 (div r ,r) 11 M TH _ 3(r) 11 u h~^(t) 

Lemma 2.3 The operators 'yrj and 7^ are linear and continuous from "^^(f^R 3 ) to 
TL 2 (r) and they can be extended to continuous linear operators from H(curl,il) and 
H (curl, O) n H (curl curl, Q), respectively, to TH" (divp, T). 

For u G Hi oc (curl, Q c ) and v G Hi oc (curl curl, f2 c )) we define 7^,11 and 7^-v in the 
same way and the same mapping properties hold true. 

Recall that we assume that the boundary T is smooth and topologically trivial. For 
a proof of the following result, we refer to [31 IHl 124] . 

Lemma 2.4 Let t£l. The Laplace- Beltrami operator 

Ar = divr Vr = — curlr curlr (2.4) 

is linear and continuous from H t+2 (T) to i/*(r). 

It is an isomorphism from i7* +2 (r)/R to the space Hl(T) defined by 

ueHt{T) ueH\T)and f u = °- 

This result is due to the surjectivity of the operators divr and curlr from TH t+1 (r) to 
Hl(T). 

We note the following equalities: 

curlr Vr = and divr curlr = (2-5) 
divr(n x j) = — curlr j and curlr (n x j) = divr j (2.6) 
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2.2 Pseudo-homogeneous kernels 

In this paper we are concerned with boundary integral operators of the form : 
/Cru(x) = vp. / k(y, x - y)u{y)do{y), x 6 T 



(2.7) 



where the integral is assumed to exist in the sense of a Cauchy principal value and 
the kernel k is weakly singular, regular with respect to the variable y £ T and quasi- 
homogeneous with respect to the variable z = x — y 6 M 3 . We recall the regularity 
properties of these operators on the Sobolev spaces H S (F), s£R available also for their 
adjoints operators: 



/Cp(n)(x) = vp. J k(x, y — x)u(y)da(y), x£T. 



(2- 



We use the class of weakly singular kernel introduced by Nedelec ( |24| p. 176). More 
details can be found in [T3l [T71IT91 l22l l35l IM] . 

Definition 2.5 The homogeneous kernel k(y, z) defined onF x (R 3 \{0}) is said of class 
—m with m > if 



sup sup 

1 1*1=1 



0M 0101 



dy a dz? 



k(y,z) 



< C a /or multi-index a and f3, 



dzi 3 



k(y, z) is homogeneous of degree — 2 with respect to the variable z 



for all |/3| = m and D™k(y,z) is odd with respect to the variable z. 

Definition 2.6 The kernel k G < ^ ° (r x (M 3 \{0})) is pseudo-homogeneous of class —m 
for an integer m such that m ^ 0, if for all integer s the kernel k admit the following 
asymptotic expansion when z tends to 0: 



N-l 



k(y, z) = k m (y, z) + ^ k m+j (y, z) + k m+N (y, z), 



(2.9) 



where for j = 0, 1, ...,N — 1 the function k m+ j is homogeneous of class —{m + j) and N 
is chosen such that k m+ N is s times differentiables. 

For the proof of the following theorem, we refer to |24j . 

Theorem 2.7 Let k be a pseudo-homogeneous kernel of class —m. The associated oper- 
ator /Cr given by (|2.7p is linear and continuous from H S (F) to H s+m (F) for all s € M. 

We have similar results for the adjoint operators /Cp. 
The following theorem is established in |13j . 



7 



Theorem 2.8 Let k be a pseudo-homogeneous kernel of class —m. The potential operator 
V defined by 

V{u){x) = [ k(y,x- y)u(y)da(y), x G R 3 \T (2.10) 



is continuous from H s 2 (r) to H s+m (Q) U H^ m (Q c ) for all positive real number s. 
2.3 The electromagnetic boundary integral operators 

We use some well known results about electromagnetic potentials. Details can be found 

in OH El El]. 

Let k be a complex number such that Im(ft) > and let 

e in\x-y\ 

G(k, \x - y\) 



4ir\x — y\ 

be the fundamental solution of the Helmholtz equation 

An + k 2 u = 0. 
The single layer potential ip K is given by : 

(ip K u)(x) = J G(k, \x - y\)u{y)da{y) x G M 3 \r, 

and its trace by 

V K u{x) = / G(k, \x - y\)u(y)da(y) x G T. 



The fundamental solution is pseudo- homogeneous of class —1 (see I24|). As conse- 
quence we have the following result : 

Lemma 2.9 Let s G E. The operators 

V K :H S -^{T)^H S+ ^{T) 

are continuous. 

We define the electric potential ^ e k generated by j G TH~2(divr,r) by 

^eJ ■= ki/j k ] + k _1 VV> k div r j 

This can be written as *&e k ] '■= k1 curl curl tp K ] because of the Helmholtz equation and 
the identity curl curl = — A + Vdiv (cf [3]). 

We define the magnetic potential V?m k generated by m G TH~2 (divr, T) by 

V^M^m := curiam. 

We denote the identity operator by I. 



Lemma 2.10 The potentials \& e k st^M K are continuous from TH 2(divr,r) to Hj oc (curl 5 
For] € TH~2(divr,r) we have 

(curl curl -k 2 1)^ e J = and (curl curl -K 2 I)^M K m = in R 3 \T 

and ^ e k ] an d ^M K m satisfy the Silver- Mutter condition. 

We define the electric and the magnetic far field operators for j € TH~2~(divr,r) and an 
element x of the unit sphere S 2 of M 3 by 



(2.11) 



*M Ke m = inxx(J e- lKiy \{y)da{y)^ . 

These operators are bounded from TH - 2~(divr,r) to T ( c tf°° (S 2 )) 3 . 
We can now define the main boundary integral operators: 

C K = -i{ 7 z> + 7d}*e« = ~\{lN + 1 c n }^m k , 

M K = -\{ 1D +7£)}*M K = -\{lN+l C N }^E K - 

These are bounded operators in TH 2 (divr, T). We have 

C K j(x) = -kJ^ n(x) x (G(k, \x - y\)](y))da(y) + k~ x J curlf (G(k, \x - y\) div r \{y))da{y) 
= [-k n x V K j + k' 1 curl r V K div r j) (x) 

and 

M K j(x) = ~y n(x) x curF(G(K, \x - y\)'](y))da(y) 
= (D K j-B K j)(x), 

B K j(x) = J V x G(k, \x - y\) (}(y) ■ n(x)) da(y), 

D K j(x) = J (V x G(k, \x - y\) • n (x))j(y)da(y). 

The kernel of D K is pseudo- homogeneous of class —1 and the operator M K has the same 
regularity as D K on TH~2 (divp, T), that is compact. 

We describe briefly the boundary integral equation method developped by the autors 
[9] to solve the dielectric scattering problem. 

Boundary integral equation method : This is based on the Stratton-Chu formula, 
the jump relations of the electromagnetic potentials and the Calderon projector's formula 
(see MM)- 



with 
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We need a variant of the operator C K defined for j £ TH 2 (divr, T) by : 

C j = n x Vq j + curlr Vb divr j. 

The operator Cg is bounded in TH 2 (divr, T). We use the following ansatz on the 
integral representation of the exterior electric field E s : 

E s = - *e Kc j - iV *u Kc Co j in R 3 \n (2.12) 

r\ is a positive real number and j £ TH ^ (divr, T). Thanks to the transmission conditions 
we have the integral representation of the interior field 

Ei = -^E Ki {Y N E inC + Nej}) - (*M Ki h c D E mc + L e j}) in SI (2.13) 

where p = and 

N e = Ql-M^ + ir, C Ke C* . 

We apply the exterior Dirichlet trace to the righthandside (|2.13p . The density j then 
solves the following boundary integral equation: 



s j = p + Mk Lcj + ° Ki Nc j = ~ p + Mk io£inc + Ck; 7jv - £mc sur r - 

The operator S is linear, bounded and invertible on TH~2 (divr, T). 

If we are concerned with the far field pattern E°° of the solution, it suffices to re- 
place the potential operators ^e k an d ^M K by the far field operators and Vl/j§ 
respectively. 

The solution E(il) = (E*($7), E s ($7)) and the far field pattern E°°(f2) consists of ap- 
plications defined by integrals on the boundary T and if the incident field is a fixed data, 
these quantities depend on the scatterrer SI only. 



3 Some remarks on shape derivatives 

We want to study the dependance of any functionals F with respect to the shape of the 
dielectric scatterer SI. The fi-dependance is highly nonlinear. The standard differential 
calculus tools need the framework of topological vector spaces which are locally convex at 
least [33j, framework we do not dispose in the case of shape functionals. An interesting 
approach consists in representing the variations of the domain SI by elements of a function 
space. We consider variations generated by transformations of the form 

+ r(x) 
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of any points x in the space M 3 , where r is a vectorial function defined (at least) in the 
neiborhood of f2. This transformation deforms the domain Vt in a domain $7 r of boundary 
T r . The functions r are assumed to be a small enough elements of a Frechet space X in 
order that (I + r) is an isomorphism from V to 

r r = (I + r)T = {x r = x + r (x) ; x G T} . 

Since we consider smooth surfaces, in the remaining of this paper, the space X will be 

the Frechet space ^(R 3 ,R 3 ) = p| ^(R 3 ,R 3 ) undowed with the set of non decreasing 

fceN 

seminorms (|| • ||fc)fceN where ^(M^R 3 ) with k G N is the space of £;-times continuously 
differentiable functions whose the derivatives are bounded and 



|r||jfc = sup sup 

0<p<fc xGR 3 



Mix) 



For e small enough we set 



=|r£ (V°°(1R 3 )) 3 , doo(0,r) < e 

where doo is the metric induced by the seminorms. 
We introduce the application 

We define the shape derivative of the functional F trough the deformation $7 — > as 
the Gateaux derivative of the application Fq in the direction £ € Af. We write: 

DJ TO<1 -!,„*>«>• 

3.1 Gateaux differentiability: elementary results 

Frechet spaces are locally convex, metrisable and complete topological vector spaces on 
which we can extend any elementary results available on Banach spaces. We recall some 
of them. We refer to the Schwarz's book f 33 1 for more details. 



Let X and 3^ be Frechet spaces and let U be a subset of X. 

Definition 3.1 (Gateaux semi-derivatives) The application f : U — > y is said to 
have Gateaux semiderivative at tq G U in the direction £ G X if the following limit exists 
and is finite 

-/Fo;£] = hm : = 7T7 1 /(ro + *£). 



a r JLU ' SJ t^o t dt 



t=0 



Definition 3.2 (Gateaux differentiability) The application f : U — > y is said to be 
Gateaux differentiable at ro G U if it has Gateaux semiderivatives in all direction £ G X 
and if the map 

is linear and continuous. 
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We say that / is continuously (or c <f 1 -) Gateaux differentiable if it is Gateaux differentiable 
at all tq G U and the application 

^-f:(r ;O£UxX^-^f[r ;t]ey 

is continuous. 

Remark 3.3 Let us come to shape junctionals. In calculus of shape variation, we usually 
consider the Gateaux derivative in r = only. This is due to the result : If J-q is Gateaux 
differentiable on then for all £ G X we have 

— ? a [r ; £] = DF(n ro ;£ o (I + r ) _1 ) = —Fn rQ [0; £ o (I + ro)" 1 ]. 

Definition 3.4 (higher order derivatives) Let m G N. We say that f is (m + 1) -times 
continuously (or c to m+1 -) Gateaux differentiable if it is c € m -Gateaux differentiable and 

d m 

is continuously Gateaux differentiable for all m-uple (£i,---,£ m ) G X m . Then for all 
ro G U the application 

(6, . . . G m- ^-^-/[roja, • • • ,£m+i] G y 

is (m + 1) -linear, symetric and continuous. We say that f is 'tf 00 -Gateaux differentiable 
if it is c € m -Gateaux differentiable for all m G N. 



We use the notation 



am am 

fM = ^ t J(ro + tO. (3.1) 



dr m 1 1 dt 
If it is ^ m -Gateaux differentiable we have 

^/[r„, 6, ■ ■ ■ , U = - { E(-!) m ^ £ & + ■■■ + &,]• (3-2) 

p=l l<ii <---<jp<m 

To determine higher order Gateaux derivatives it is more easy to use this equality. 

The chain and product rules and the Taylor expansion with integral remainder are still 
available for ^""-Gateaux differentiable maps ( [33j p. 30). We use the following lemma 
to study the Gateaux differentiability of any applications mapping r on the inverse of an 
element in a unitary topological algebra. 

Lemma 3.5 Let X be a Frechet space and y be a unitary Frechet algebra. Let U be an 
open set of X. Assume that the application f : U — > y is Gateaux differentiable atrQ^iU 
and that f(r) is invertible in y for all r G U and that the application g : r h-> f(r)~ 1 
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is continuous at tq. Then g is Gateaux differentiable at r$ and its first derivative in the 
direction £ G X is 

^fM = -firor 1 o |;/[ro,£] o /(ro)" 1 . (3.3) 
Moreover if f is C € m -Gateaux differentiable then g is too. 

PROOF. Let £ G <Y and t > small enough such that (r + i£) G 17, on a: 

g(r + tO - g(r ) = /(ro)" 1 o /(r ) o /(r + *£)-* - /(ro)" 1 o /(r + i£) o /(r + *£)-* 
= /(ro)" 1 o (/(r ) - /(r + i£)) ° /(r + *0 _1 
= /(ro)" 1 o (/(r ) - /(r + i£)) ° /(ro)" 1 

+/(ro)" 1 o (/(r ) - /(r + i£)) ° (/(r + *£)-* " /(ro)" 1 ) • 

Since g is continuous in ro, we have lim (/(ro + 1^)" 1 — /(ro) -1 ) =0 and since / is 
Gateaux differentiable in ro we have 

/(ro)- 1 o(/(ro)-/(ro + ^))o/(r )- 1 x 8 , 

Jim = -(/M) o ^/[r ,£] o (/(r )) . 

As a consequence 

,. g(ro + tO - g(r ) i d x 

t™ t = -(/(ro)) ^/Fo,C]°(/(ro)) • 



4 Gateaux differentiability of pseudo-homogeneous kernels 

Let x r denote an element of T r and let n r be the outer unit normal vector to T r . When 
r = we write no = n. We note again da the area element on T r . 

In this section we want to study the Gateaux differentiability of the application map- 
ping r G to the integral operator /Cr r defined for a function u r G H s (T r ) by: 

K Tr u r {x r ) = vp. / k r (y r ,x r - y r )u r (y r )da(y r ), x r G T r (4.1) 

and of the application mapping r G to the potential operator V r defined for a function 
u r G H s (T r ) by: 

V r u r (x) = k r (y r ,x - y r )u r (y r )da(y r ), x G K, (4.2) 
where k r G (r r x (M 3 \{0})) is a pseudo-homogeneous kernel of class — m with m G N. 
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We want to differentiate applications of the form r i— > J~u(r) where the domain of 
definition of J-fi(r) varies with r. How do we do? A first idea, quite classical (see 
|25 [ 127 1 1^9])- is that instead of studying the application 

r € fl e °° Jh(r) E ^ k {T r ) 

we consider the application 

r € 5,°° ^ 7h(r) o (I + r) E ^ fc (r). 

An example is r i— >• n r . This point of view can be extended to Sobolev spaces H S (T), 
s£l. From now we use the transformation r r which maps a function u r defined on T r 
to the function u r o (I + r) defined on T. For all r E -B^°, this transformation r r admit an 
inverse. We have 

(r r u r )(x) = u r (x + r(x)) and (t~ u)(x r ) = u(x). 
Then, instead of studying the application 

r e s e °° _> /c r ,. e if c (tf s (r r ), #*+ m (r r )) 

we consider the conjugate application 

r6 5 e °°4 r./C^r- 1 E if c (iT s (r), #* +m (r)) . 

In the framework of boundary integral equations, this approach is sufficient to obtain the 
shape different ability of any solution to scalar boundary value problems \27\ I29| . 
Using the change of variable x i— > x r = x + r(x), we have for u E H S (T): 

^ W -X^ + ^- + K.,-.-^)*^).. 6 r 

where J r is the jacobian (the determinant of the Jacobian matrix) of the change of variable 
mapping x E T to x + r(x) E T r . The differentiablility analysis of these operators begins 
with the jacobian one. We have 

J r = Jacr(I + r) = |[u> r [| with oj r = com(I + Dr| r )no = det(I + Dr| r ) T (I + Dr| r ) _1 n, 
and the normal vector n r is given by 



r- 1 



Co',. 



The first derivative at r = of these applications are well known |12[ [25] . Here we present 
one method to obtain higher order derivative. 



Lemma 4.1 The application J mapping r E to the jacobian J r E c rf QO (T,M) is 
Gateaux differ •entiable and its first derivative at ro is defined for £ E ^"(R^R 3 ) by: 

d,J 

frfa't] = J ro(i~r divr rQ T~)£. 
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Proof. We just have to prove the ^^-Gateaux differentiability of W : r \— > w r . We 
do the proof for hypersurfaces T of W n , n £ N, n > 2. We use local coordinate system. 
Assume that T is parametrised by an atlas (Oj, (f>i)i<i< p then T r can be parametrised by 
the atlas (I + r) o (f>i)\<i< p . For any iGT, let us note ei(x), e2(x), . . . , e n _i(x) the 
vector basis of the tangent plane to V at x. The vector basis of the tangent plane to T r 
at x + r(x) are given by 

ej(r, x) = [(I + Dr)(x)]ei(x) for i = 1, . . . , n — 1. 



n-l 



A ei(r,x) 



Thus, we have u) r {x) 



i=i 



Since the applications r i— > ej(r, x), for i = 1, . . . , n— 1 



n-l 

A ei(x) 
i=l 

are < £?°°- Gateaux differentiable, the application W is too. Now want to compute the 

and t small enough. We have at 



derivatives using the formula (|3.2p . Let £ £ ^ 

r G B, 



oof wan Trpn\ 



n-l 

A (I + Dr + tD£)ei(x) 

i=l 



&t r ' 



t=0 



n-l 

A ei(a!) 

i=l 



To simplify this expression one have to note that 



[D£0)h(aO = [De(x)][(I + Dr )(x)]- 1 [(I + Dr )(x)]e 4 (x) 

= [D£(x)][D(I + ro)" 1 ^ + r (x))][(I + Dr )(x)] ei (x) 
= [( T r D^ 1 )C(^)]e t (r ,a;) = [(r ro Dr^r" 1 ^ (z)]ei(r , x) 

NB: given a (n x n) matrix j4 we have 



n— 1 



n-l 



• • • x ej_i x j4ej x e$+i x • • • = (Trace(j4)I — T A) A ej. 



i=l 



i=l 



Thus we have with A = [r ro D r r" 1 ^ and 5 = I, #i(A) = Trace(A)I - T A 

= «^ro( r ro rl T- )) 



(#) 



dW 



Or 
d m W 



r o,C] = Jr ((7"r divr rQ r^ 1 )^ ■ r ro n ro 
= [^(AJdWCro), 

r ,£] = [B m (A)Z]W(r ) 

= E("l) m 7 ^[BiCA^B^CAJ^WCro) for m = 1, . . . ,n - 1, 

i=i (m-t)! 



I <9r m 



[ r 0j£] = for all m>n. 
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It follows that 

t^ r °' el = \\w\r )f W [r °^ ] ' W{ro) = ^ [r °^ ] ' Tro " ro = Jro(rro diVrv ° Tr « 1)e 



Thanks to (#) we deduce easily the Gateaux differentiability of r i-> r r n r . 

Lemma 4.2 TTie application N mapping r € -B^° to r r n r = n r o (I + r) G ^°°(r,lR 3 ) is 
< ^ 700 Gateaux- differentiable and its first derivative at ro is defined for £ € ^°°(1R 3 , M 3 ) 6y: 

^h,C] = -[r ro Vr, o r- 1 £]AT(ro). 

PROOF. This results from the precedent proof and we have 

dN 1 dW 1 f dW \ 

= - [^oVr^r-^] 7V n ro . 



To obtain higher order shape derivatives of these applications one can use the equalities 
(#)and 

1 1 Tf n f 1 1 1 <j 

(*) { d m N-N 

, dr m M = Ofor allm> 1. 

As for exemple we have at r = in the direction £ G ^^(M 3 , M 3 ): 

^[0,e]=div r eand^[0,£] = -[V r e]n, 

d 2 J 

^2 [0,6,6] = -TVace([V r 6][Vr6]) + div r 6 • div r 6 + ([V r 6]n ■ [V r &]n) • 
Notice that Trace([V r 6] [Vr6]) = Trace([V r 6][V r 6D, 
d 2 N 



dr 2 



[0,6,6] = [Vr6][Vr6]n + [V r 6][Vr6]n - ([V r 6]n ■ [V r 6]n)n. 



For n > 3 it needs too long calculations to simplify the expression of the derivatives and 
we only obtain the quadratic expression. In the last section we give a second method to 
obtain higher order derivatives using the Gateaux derivatives of the surface differential 
operators. 
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Remark 4.3 We do not need more than the first derivative of the deformations £. As a 
consequence for hyper surf aces of class c £ k+1 , it suffice to consider deformations of class 
to conserve the regularity of the jacobian and of the normal vector by differentiation. 



The following theorem establish sufficient conditions for the Gateaux differentiabil- 
ity of the boundary integral operators described here above and that we obtain their 
derivatives by deriving their kernels. 

Theorem 4.4 Let k G N. We set (T x T)* = {(x,y) G T x T; x / y}. Assume that 

1) For all fixed (x,y) G (T x T)* the function 

f : B™ -»• C 

r h-> k r (y + r(y),x + r(x) - y - r(y))J r (y) 

is c f k+1 -Gateaux differentiable. 

2) The functions (y, x — y) t- > f(ro)(y, x — y) and 

(y,x-y) ^ - d -^f[ro,Zi,...,£,a]{y,x -y) 
are pseudo-homogeneous of class —m for all ro G B%° , for all a = 1, . . . , k + 1 and for all 

6, •••,6+1 e^T r::i Tp;; 



Then the application 



B™ -> ^ c (H s (T),H s+m (T)) 

r h+ T r K,T r T~ l 

is ^ -Gateaux differentiable and 

d k f d k 

— f {T r /Cr r T r -1 } [r , 6, • • • , t,k]u{x) = ^/fa), £i, •••,£k](v,x- y)u{y)da(y). 

PROOF. We use the linearity of the integral and Taylor expansion with integral remain- 
der. We do the proof k = 1 only. Let ro G B^°, £ <G ^°°(]R n , R n ) and t small enough such 
that r + i£ G . We have 

f(ro+t£, x, y)-f(r , y, x-y) = t— [r , f](y, x-y)+t 2 (l-A)^- [r +Ai£, f ] (y, x-y)d\. 

We have to verify that each terms in the equality here above is a kernel of an operator 
mapping H S (T) onto .fP +m (T). The two first terms in the left hand side are kernels of 

d 2 f 

class — m and by hypothesis q~~2"[0,£] * s a ^ so a kernel of class — m. It remains to prove 
that the operator with kernel 

f 1 d 2 f 
(x,y)^ J (l-X)^±[ ro + \tC,C](x,y) 
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d 2 f 

acts from H S (T) to H s+m (T). Since tto^o + is pseudo-homogeneous of class — m 

or z 

for all A G [0, 1], it suffice to use Lebesgue's theorem in order to invert the integration 
with respect to the variable A and the integration with respect to y on T. 



jf (l-A) A/i, .(](,('.;/)// (y)(/(T(y) j r/A 



H s + m (r) 



< s u PAe[o,i] 

< C\\u\\ H s (r) . 
We then have 



dr 2 

d 2 f 



r dr 2 



[ro + \t£,$,](x,y)u(y)da(y) 



ifs+m(r) 



j + ^' x ' y) u (y) da (y) ~ J r f( r o> x > y) u (y) d(j (y) 



dj_ 

dr 



[ro,£,](x,y)u(y)da(y) + t 



my 



d 2 f 

A )^i [ r ° + xt a^,]{x,y)d\ ) u (y)d°(y)- 



We pass to the limit in t = and we obtain the first Gateaux derivative. For higher order 

Qk 

derivative it suffice to write the proof with Tr-r/ko; £i> ■ ■ ■ , £fc] instead of /. The linearity, 

or K 

the symetry and the continuity of the first derivative is deduced from the kernel one. ■ 

Now we will consider some particular classes of pseudo-homogeneous kernels. 
Corollary 4.5 Assume that the kernels k r are of the form 

kriyr-i %r Vr) = G{x r y r ) 

where G is pseudo-homogeneous kernel which do not depend on r. Then the application 



-> ^ c {H s {T),H s+m {T)) 



is < ^°° -Gateaux differentiable and the kernel of the first derivative at r = is defined for 
£ G ^(R 3 ,R 3 ) by 



d {G(x + r(x) - y - r(y))} 
dr 



[o, ei = (a^) - ay)) ■ v z g(* - y ) + g{ X - y ) di Vr 



PROOF. For all fixed (x,y) G (r x T)* , consider the application 

/ :U i-> /(r, x, y) = G(a: + r(x) - y - r(y))J r (y) G C. 

We have to prove that r i— >■ f(r) is ^^-Gateaux differentiable and that each derivative 
define a pseudo-homogeneous kernel of class — m. 
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> Step 1: 

First of all we prove that for (x,y) S (r x T)* fixed the application r i— >• f(r,x,y) is 
infinitely Gateaux differentiable on . By lemma 14.11 the application r i— > J r (y) is 
infinitely Gateaux differentiable on B^°, the application rHi + r(x) is also infinitely 
Gateaux differentiable on 2^° and the kernel G is of class < ^ , °° on IR 3 \{0}. Being composed 
of applications infinitely Gateaux differentiable, the application r i— > f(r,x,y) is too and 
using Leibniz formula we have : 

Qk 

Q-jf[ro,€i,---,€k](x,y) = 

k Qa d k ~ a J r (y) 

E E {G{x + r(x) -y- r(y))} [r ,^ (1) , . . .,g CT ( g )] fc _ r a [r ,^( a+ i),. . • ,^(fc)] 

where S^T denote the non decreasing permutations of {1, ... , k} and 

Qa 

{G{x r - y r )} [r ; £i, . . . ,£ Q ] = D%G[x+r (x)-y-r (y);£i(x)-£i(y), . . . , € a (x)-£ a (y)]. 



dr a 

> Step 2: 

We then prove that each derivative define a new pseudo-homogeneous kernel of class — m 
that is to say that for all k G N and for all A;-uple . . . the application 

Qk 

d k ~ a J r 

is pseudo-homogeneous of class — m. Since — k _ a [ro,£i, ■ ■ ■ £ ^ > °°(r,R) we have 

to prove that 

^ ^ { G ( x + r(x ) ~ y ~ r (y))l [ r o,6,- • • ,&*] 

defines a pseudo-homogeneous kernel of class —m. By definition, G(z) admit the following 
asymptotic expansion when z tends to zero: 

N-l 

G{z) = G m {z)+ Y,Gm +j (z) + G m+N(h, z) 

3=1 

where G m +j is homogeneous of class — (m + j) for j = 0, . . . , N — 1 and G m +7\r is of 
arbitrary regularity. Using Taylor formula we obtain the following result : 

Proposition 4.6 Let G m (z) be an homogeneous kernel of class —m and any deformations 
£ = (&,...,&»)€ ^°°(M 3 , R 3 ). TTie function 

(x,y) ^ D a G m [x + r (x) - y- r (y);£i(x) - . . -,£ a {x) - £ a (y)] 

is pseudo-homogeneous of class —m. 
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The application mapping (£1, . . . , £ a ) £ , R ) to the integral operator of kernel 

d k {G(x + r{x) -y-r(y))}. 



is clearly linear and continuous for all rg € B^ 



■[r ,£i, . ..,6] 



Example 4.7 (Single layer kernel) We note V£ the integral operator defined for u r £ 
H s (T r ) by 

V*u r (x) = / G(k, \x r - y r \)u r (y r )da(y r ). 



rv 

The application 

b s -> J£ c (h s (r) , h s+1 (r) ) 

zs ff 00 -Gateaux differentiate and its first derivative atr = in the direction £ € < ^ >t 
is 

rr_— 1 



where in M 3 «;e /iat>e 



' jr ' 1 ' 7 [0, £]«(*) = / k'(y,x- y)u{y)da{y) (4.3) 



r 



(x, y) =G(k, \x - y\) . , in - , + div r £{y) 

V \x-y\ V \ x -y\J 

Example 4.8 (Double layer kernel) We note D r K the integral operator defined for u r G 
H s (T r ) by 

D r K u r (x) = / n r (x r ) -V z G(K,\x r - y r \)u r (y r )da(y r ). 

The application 

B s -> ^ c {H*{T),H s + l {T)) 



r ^ T r D r K r; 



is Gateaux- differentiate . 
Proof. We have 

n r (x r ) ■ V z G(k, \x r - y r \)u r (y r ) = n r (x r ) ■ (x r - y r ) G ^ — f~~ (in - -, — ^ — 

\x r y r \ \ \x r y r 

We have to prove that 

r G B°° h-> ( T r"r)(a:) ; (a + r(x) - y- r{y)) 
\x + r(x) — y — r(y)\ 3 

is ^ ,0 ° Gateaux differentiable and that the derivatives are pseudo-homogeneous of class 
— 1. To do so we use local coordinates as Potthast did in |29] and prove that 

d k (r r ti r )(x) • (x + r{x) - y - r{y)) 
lro,£i,...,&J 

behaves as \x — y\ 2 when x — y tends to zero. ■ 
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Each domain Q, is a countable union of compact subset of Q: f2 = \J K p . Instead of 

P >1 

studying the application 

reBfH T n (r) G if c (H s (T r ), H s+m (fl r )) 
we consider the application 

r G Bf M- JhWr" 1 € J^ c (iT (r), H s+m (K p )) . 
We use this approach for potential operators. We have for u G Lf s ~2(T) 

T> r T~ l {u)(x) = J k r (y + r(y),x - y - r(y))u(y)J r (y)da(y), x e K p . 

Theorem 4.9 Lei s G R. Lei G(-z) be a pseudo-homogeneous kernel of class — (m + 1) 
m G N. ylssume that for all r G , we have k r (y r , x — y r ) = G(x — y r ). Then the 
application 

£ e °° -> ^ c (F s -5(r),^°°(^ p )) 



r i-t> P r r r 



-i 



is infinitely Gateaux differentiable and 
d k V r T~ 1 f d k 

— Q r k [ro,ti, • • • , &]«(a0 = J r Q^k {°( x ~ V ~ r (y))^r(y)} [ r o, &,•••> CkMl/)*r(l/)- 

/is /irsi derivative at r = in the direction £ G ^"(R 3 , R 3 ) is i/te integral operator 
denoted by V 1 with the kernel 

{ G ( x - y - r ( y ))} [ ro ,£] = -fa) . v z G(x -y) + G(x - y) div r (4.4) 

T/ie operator can be extended in a linear and continuous integral operator from 

H s -h(T) to H s+m (n) and H°+ m (fl). 

PROOF. The kernel and its higher order derivatives are of class on K p . 
Since $7 is an increasing union of compact manifolds we can define a shape derivative on 
the whole domain Q. Let us look at the first derivative : the term G(x — y) divr has 
the same regularity than G(x — y) when x — y tends to zero wheareas £(y) • VG(x — y) 
loose one order of regularity. As a consequence the kernel must be of class — m + 1 in 
order that its first derivative acts from H S -^{T) to H s+m {£l) and jyf+ m (fi c ). ■ 

Remark 4.10 We conclude that the boundary integral operators are smooth with respect 
to the domain whereas the potential operators loose one order of regularity at each deriva- 
tion. We point out that we do not need more than the first derivative of the deformations 
£ to compute the Gateaux derivatives of these integral operators. 
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Example 4.11 (Single layer potential) We denote by Uj r K the single layer potential 
defined for u r G H s (T r ) by 



ip r K u r (x) = / G(k, \x - y r \)u r (y r )da(y r ), x G R \T r . 
Jr r 

The application 

r h-s> T r u) r K T- 1 

is infinitely Gateaux differentiable. Its first derivative at r = can be extended in a linear 
and continuous operator from if s ~2(r) to H s (£l) U Hf oc {VL c ). 

Since the potential operators are infinitely Gateaux differentiable far from the bound- 
ary we have the following result by inverting the derivation with respect to r and the 
passage to the limit \x\ — > oo. 

Example 4.12 Let s G R. We denote by rh^' r the far field operator associated to the 
single layer potential defined for u r E H s (T r ) by 



^ r u r (x) = / e- lKX ^u r (y r )da(y r ), x G S 2 . 

The application 

lT ,oo,r —1 

r M- W K r r 

is infinitely Gateaux differentiable and its first deriavtive atr = 0is defined foru G H S (T) 
by: 

d ^Q r T ^ [o. = ( div r av) - i*& • av)) »(i/)<fr(i/)) • 

5 Shape differentiability of the solution 

Let E mc be an incident electric field which is a fixed data. The aim of this section is 
to study the shape differentiation properties of the application E mapping the bounded 
scatterer Q to the solution E(fl) = (E l (fi), E S (S7)) G H /oc (curl, K 3 ) to the dielectric 
scattering problem by the obstacle Q lit by the incident field E mc established in section 
[TJ To do so we use the integral representation of the solution. 

We set S\r) = E*(SV) and S s (r) = E s (ft r ) and we denote C£ r , C£ and 

M£ the potential operators and the boundary integral operators on T r and 7^, 7]^- , 7^ 
the 7^ trace mappings on T r . We have : 

^tot( r j = E <nc + ^ s (r) (5.1) 

with 

<f s ( r ) = - %rfSi r M C * r ) j r dans ^ = W^^T r (5.2) 
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where j r solves the integral equation 

Sr\ „ I T i A/fT i „T c*"c /nrr _,r 

Jr = "Pi -2 I + M «i J TD E _C VYjV Ke E , 

and 

r(r) = --^.7^1 - *k.7S r ^°*(r) dans ft r (5.3) 

Recall that the operator S r is composed of the operators CT , MI , CI! et Ml and that 
these last ones are defined on the space TH~2 (divr r , I"V). 



5.1 Variations of Helmholtz decomposition 

We have to turn out many difficulties. On one hand, to be able to construct shape 

derivatives of the solution it is necessary to prove that the derivatives are defined on the 

same spaces than the boundary integral operators themselves, that is TH _ 2(divr,r) (if 

we derive at r = 0). On the other hand, the very definition of the differentiability of 

operators defined on TH _ 2(divr,r) raises non-trivial questions. The first one is : How 

to derive applications defined on the variable space TH~5 (divr r , r r )? 

A first idea is to insert the identity r~ x r r = I i between each operator in the integral 

3 r i-r2(r r ) y b 

representation of the solution in order to consider integral operators on the fixed boundary 
r only and to study the differentiability of the applications 

r i y T r C r K T T , 
r i— > t t M t k t^ x , 

but many difficulties persist as Potthast pointed out [28]. The electromagnetic boundary 
integral operators are defined and bounded on tangential functions to T r .The restriction of 
the operator t t M t k t~ x to tangential densities to r r , has the same regularity of the double 
layer potential operator. If we differentiate r r M^r~ , we will not obtain an operator with 
the same regularity than M K and acting on TH _ 2(divr,r) since: 

r r (TH-^(divr r ,r r )) + ThH(div r , T). 
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The incident field E mc is analytic in the neighborhood of T thus 7£,E mc G TH 2 (divr r , T r ) 
for all r G Bf. Set f(r) = t t (Yd^) ■ For £ G ^{M. 3 ,R 3 ), the Gateaux semiderivative 



q p Q j~ (i^^) 

— is not tangent to T anymore it follows that M K — is not defined. We 

at |t=o at |t=o 

have the same difficulties we the Neuman trace 7^ and the other operators. 

As an alternative, the idea of R. Potthast was to introduce projectors on the tangent 
planes of the surfaces V and T r . Let us note 7r(r) the orthogonal projection of any 
functions defined on T r onto the tangent plane to T. This is a linear and continuous 
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operator from the continous vector function space on T r to the the space of continuous 
tangential function to T and for u r G (^(r,.)) 3 we have 

(7r(r)u r )(x) = u r (x + r(x)) — (n(x) • u r (x + r(x))) n(x). 

Proposition 5.1 The restriction of ir(r) to the continuous and tangential functions to 
T r admit an inverse, denoted by 7r _1 (r). The application tt (r) is defined for a tangential 
function u to T by 

f -i, , w . ( f , , ,n r (x + r(x))-u(x) 

(tt (r)u)(x + r(x)) = u(x) — n(x)- 



n r (x + r(x)) ■ n(x) ' 
And we have 7r _1 (r)u G TH s (r r ) if and only if u G TH S (T). 

In the framework of the space of tangential continuous functions it suffices to insert the 
product 7r _1 (r)7r(r) = I-pg'Ofrv) i J1 the integral representation of the solution to lead us to 
study boundary integral operators defined on T^°(r) which do not depend on r anymore 
but here we would obtain operators defined on 

7r(r) (TH-5(divr r ,r r )) = {u G TH-^diviY^^u) G fH(r r )} . 



This space depends again on the variable r and do not correspond to TH~2 (divr, T). Our 
approach consist in using the Helmholtz decomposition of the spaces TH _ 2 (divr r , T r ) for 
r G and to introduce a new invertible operator P r defined on TH~2 (divp,., r r ) and 
which is not a projection operator. 

We have the following decomposition. We refer to [llj for the proof. 

Theorem 5.2 The Hilbert space TH~2(divr,r) admit the following Helmholtz decompo- 
sition: 

ThT5(div r ,r) = V r (H§(T)/R) 0curl r (h?(T)/R\ . (5.4) 

Since the real e is chosen such that for all r G the surfaces T r are still regular and 
simply connected, then the spaces TH~2 (divr r ,r r ) admit the same decomposition. 
Let j r G TH _ 2(divr r ,r r ) and let Vr r p r + curlr,, q r its Helmholtz decomposition. Since 

3 1 

p r G Hz(T r ) and q r G H^(T r ), their change of variables from T r to T, T r (p r ) and T r (q r ), 
are in H 2 (r) and H 2 (r) respectively. The following operator : 

P r : TH-5(divr r ,r r ) — > ThW (div r , T) (55) 

J*r = v r r Pr + rotr r q r >-» j = V r r r p r + rot r r r gv 

is well-defined. 

The operator P r transforms a tangential vector field j r to T r in a tangential vector 
field j to r. This operator is linear, continuous and admit an inverse P" 1 given by : 



P- 1 : TH-2-(div r ,r) — > ThH (div r „, T r ) 

j = V r p + rot r q >->• j r = V r ,, r" 1 (p) + rot r ,, r" 1 (?) 
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Obviously we have when r = that P r = P r 1 = I _ 1 We insert the identity 

TH ^ (divp ,r) 

I 1 = P~ P~ between each operator in the integral representation of the solu- 

TH-2(divr r ,rV) r ^ b ^ 

tion (& l (r), S s {r)). Finaly we have to study the Gateaux differentiality properties of the 
following applications : 



-> Jg? c (TH-3(divr,r),H(curl,iirp)) : r^W^p- 1 

Bf -> Jgf c (TH-3(divr,r),H(curl,iiC p )) 

B^ -> ^ c (TH s (div r ,r),TH-5(div r ,r)) 

B f ^ c (TH-5(div r ,r),TH-5(div r ,r)) 



r^M/^P," 1 

r !-»• Pr-M^P" 1 

r ^ P^P" 1 



(5.7) 



where if„ is a compact subset of M 3 \r. 



Now let us look at the integral representation of these operators . 
> Integral representation of W E P~ x . 

The operator ^ r E P r 1 is defined for j = Vr p + curlp q € TH~2 (divp, T) and x € K p 
by: 

^P- 1 .^) = k G(k, \x - y r \) (V^r" 1 ?) (y r )da(y r ) 

+ n l G(k, \x - y r \) (rotr r r r _1 g) (y r )da(y r ) 

JY r 

- n~ l V / G(k, \x - y r \) (Ar.r^p) (y r )da(y r ). 

JY r 

D> Integral representation of W M P~ . 

The operator f ^P r _1 is defined for j = Vr p + curlp q £ TH~2(divr,r) and x £ K p 

by: 

*k p r = curl / \ x ~ Vr\) {Vr r T~ l p) (y r )da(y r ) 

+ curl / G(k, \x - y r \) (curlr r r~ l q) (y r ){y r )da{y r ). 

JY r 

D> Integral representation of P r C£P~ . 

Recall that for j r € TH" (divr r , IV), the operator C r K is defined by 

C r J r (x T ) = -Kn r (x r ) x / G(k,\x t - y r \)\ r (y r )da(y r ) 

JY r 

-/t _1 n r (x r ) x / G(k, \x r - y r \) div Tr i r (y r )d<j(y r ). 

JY r 

We want to write j r of the form Vr r -P r + curlrv Q r - Using the formula (|2.5D - (|2,6p we 
deduce that : 

divr r j r = A Fr P r et curlrv C T K j r = -Ar r Q r - 



25 



As a consequence we have for x r G T r : 

Pr(x r ) = -k Ap^ divr r ^n r (x r ) x G(k, \x r - yr\)\ r {yr)da{y r )^ (5.8) 



and 



Qr(x r ) = -k (-Ap^curl^ ^n r (av) x G(k, \x r - yr\)j r (yr)da(y r )j 

-K -1 (-Ar I .)curlr T .(-curlr r ) / G(k, \x r - y r \) div Tr } r (yr)da(y r ), 
= K A r! curl rv (^r(xr) x J G(k, \x r - y r \)\ r {yr)d(T{y r )^ 

+k~ 1 \ G(k, \x r - y r |)divr r } r (y r )da(y r ). 

The operator P r C r K P~ l is defined for j = Vr p + curlp q € TH™ 2 (divr, T) by: 

PrClP; 1 = V r P(r) + curlr Q(r), 

with 

P(r)(x) = -k (r r A^ div Fr r' 1 ) ^(r r n r )(x) x r r G(k, \ ■ -y r \)(^r r r~ 1 p)(y r )da(y r ) 

+ J r I ' -yr\)( rot r r T r 1( l)(yr)d(j(y r )^ (x) 

and 

Q(r)(x) = k (r r A^ curlr r r" 1 ) ^(r r ti r )(x) x r r G(k, | • -yr|)(Vr r r~ V)(yr)dcr(y r ) 

+^ G(k, I • -yr\)(rot rr T r r 1 q)(y r )da(y r )^ (x) 

+K~ 1 r r ^ G(k, I • -y r |)(Ar r T~V)(y r )dcr(?/ r )^ (x). 

t> Integral representation of Pj-M^P^ 1 . 
Recall that for all j r G TH~5 (divr r , T r ), the operator M£ is defined by 

M r J r (x r ) = / ((X7 Xr G(K, \x r - y r \)) ■ n r (x r ))\ r (y r )da(y r ) 
- I V Xr G{K, \x r -y r \)(n r (x r ) ■ } r (y r )) da(y r ). 

JT r 
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Using the equalities (|2.6p and the identity curl curl = —A + Vdiv, we have 

divr r M r K \ r {x r ) = n r (x r ) • / curlcurF r (G(k, \x r - y r \)]r{y r )) da(y r ) 
= K 2 n r (x r ) ■ / (G(K,\x r -y r \)j r (y r ))da(y r ) 
d f 

+ TT~ / (G(K,\x r -yr\)drvr r jr(yr))da(y r ) 

Proceeding by the same way than with the operator P r C£P^~ , we obtain that the operator 
PrMJIP- 1 is defined for j = V r p + curl r q € TH _ 5 (div r , T) by: 

PrM^p- 1 ] = VrP'(r) + rotrQ'(r), 



with 



P'{r){x) = (rrAfV- 1 ) r r |k 2 ^ n r ( • ) • {G(k, \ ■ -y r \) curlr r T^ l q{y r )} da{y r ) 

+ k 2 n r ( • ) • I • -yr|)Vr r r r _1 p(y r )} dcr(y r ) 

•/rv 



+ j£ ^ G(k, I • -y r j)(A r ,.r r 1 p)(y r )da(y r ) \ (x), 



and 



= (r r A r * curlr r r r x ) r r jjf (VG(«, | • -y r |) • n r ( • )) (rotr r r r l q){y r )da{y r ) 
+ / ((VG(k, I • -y r |) • n r ( • )) {Vv r r- X p){yr)da{Vr) 
- / VG(k, [ • -y r .|) (n r ( • ) • (rotryT" 1 ?)^)) da(y r ) 

- VG(«, [ • -y r |) (n r ( • ) • (Vr I .r- 1 p)(yr)) <M</r) j (x). 

These operators are composed of boundary integral operators with weakly singular and 
pseudo-homogeneous kernels of class -1 and of the surface differential operators defined 
in section [2j By a change of variables in the integral, we then have to study the differen- 
tiability properties of the applications 



r i— > T r Vr r T r 1 

t i y T r curlr r t" 1 

r H> T r divr r t" 1 

r h> T r curlr r r" 1 
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5.2 Gateaux differentiability of the surface differential operators 

Lemma 5.3 The application 

G : ->• if c (iF +1 (r), H s (r)) 

r i-> 7V.Vry.T~ 1 

is < tg" x> -Gateaux differentiable and its first derivative at ro is defined for £ € ^°°(IR 3 ,]R 3 ) 
— [r ,C]n = -[G(r )e]G(r ) U + (G(r )n- [G(r )C]iV(ro))iV(ro). 

Remark 5.4 iVofe i/iai we can write [ro,£] = — [G(ro)QN(ro). Since the first deriva- 
tive of N and G can be expressed in function of N and G we obtain the Gateaux derivative 
of all order iteratively. 

PROOF. In accordance to the definition (|2,ip and the lemma 14,21 to prove the c € ca - 
Gateaux differentiability of G we have to prove the < ^°°- Gateaux differentiability of the 
application 

/ : r G Bf ^ L ^ r T (WC) ^ | G J? C (H S+1 (T), H S (T)). 
Let x G r, we have 

r r ( Vr r _1 u J (x) = V (u o (I + r) _1 )| r (sc + r(x)) = T (I + Dr)^ 1 (x + r(x)) o Vu| r (x), 
and 

(I + DrV 1 (x + r(x)) = [(I + Dr) |r (x)] _1 . 



The application g : r G -Bf 3 >->■ (1+ Dr)i r G ^""(r) is continuous, and ^°°-Gateaux 

vanish. One can show that the application h:r G -B^° i— >• {a; i— )■ [g , (r)(x)] _1 } G ^^(r) is 
also Gateaux-differentiable and that we have at rg and in the direction £: 



differentiable. Its first derivative is — <?[(), £] = [D£]| r and its higher order derivatives 



~~"[ r °>£] = -h(r ) o Q^[ro,£] ° h(r ) = -h(r ) o [Df]| r o h(r ). 



and 



g^[ro,Ci, ■ ■ ■ ,&J = ^ (I + Dro)- 1 o [t^Dt- 1 ^!)] o . . . o [^Dr" 1 ^] 

where is the permutation groupe of {1, . . . , n}. Finally we obtain the < ^°°— Gateaux 
differentiability of / and we have 

df 

-QpTo^\ u = -[f( r o)C]f(ro)u. 
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To obtain the expression of the first derivative of G we have to derive the following 
expression: 

G{r)u = (r r Vr,.T r r 1 n) = r r V yr^u) — ^T r n r • ^r r V y^^^jj) Trl1r 
= f(r )u-(f(ro)u-N(r ))N(r ). 
By lemma H~2l and the chain and product rules we have 
BC 

-frlro,® = -[f(roW(ro)u + ([f(roW(ro)u-N(r ))N(r ) 

+ (f(r )u • [G(r )£]iV(ro)) N(r ) + (f(r )u ■ N(r )) [G(r )£]N(r ) 

We had the first two terms in the right handside, it gives : 
BC 

^rh),£] = -[G(r )£]f(r )u + (f(r )u ■ N(r )) [G(r )£]N(r ) 
+ (f(r )u-[G(r )Z}N(r ))N(r ) 
= -[G(r )£]G(r Q )u + (f(r )u ■ [G(r )€]N(r )) N(r ). 

To conclude it suffice to note that (f(r )u ■ [G(r )g\N(r Q )) = {G(r )u ■ [G(r )£]N(r Q )). 
■ 

Lemma 5.5 The application 

D : B™ ^f c (H s+1 (T,R 3 ),H s (T)) 
r H> r r divr r r^ 1 

is 'to 00 -Gateaux differentiable and its first derivative at ro is defined for £ € < ^j 5 °(R 3 ,R ) 
by 

^[r ,£]u = -Trace([G(r )C][G(r )u]) + ([G(r )u]tf(r ) ■ [G(r )£]iV(ro)) . 

PROOF. For u 6 H s+1 (r,M n ) we have D(r)u = Trace(G(r)u). Then we use the differen- 
tiation rules. 

Remark 5.6 Since the first derivative of D is composed of G and N and the first deriva- 
tive of J is composed of J and D, we can obtain an expression of higher order derivatives 
of the jacobian iteratively. 



Lemma 5.7 The application 

R: Bf -> j^ c (# s+1 (r),H s (r)) 

r i V r r curlrv r^ 1 

is ff 00 -Gateaux differentiable and its first derivative at ro is defined for £ € ^"(IR 3 , IR 3 ) 
by 

— [r ,Z]u= T [G(r )£]R(r )u - D(r )£ ■ K(r )u. 
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PROOF. Let u € H S+1 (T). By definition, we have R(r )u = G(r )u x N(r ). By 
lemmas 14.21 and 15.31 this application is c to°° Gateaux differentiable. We have in tq and in 
the direction £ E c tf°°(r, 



— [r ,£]u = - T [G(r )Z]G(r )u x N(r Q ) - G(r )u x [G(r )Z]N(r ). 

NB: recall that given a (3 x 3) matrix A and vectors b and c we have 

Abx c + bx Ac = Trace(^)(6 x c) - T A(6 x c). 

We deduce the expression of the first derivatives with A = — [G(ro)£], b = G(tq)u et 
c = JV(r ). ■ 

Lemma 5.8 The application 

R: Bf -> ^ c (H s+1 (r),^ s (r)) 

r h-> r r curlry r" 1 
is 'lo 00 -Gateaux differentiable and its first derivative at ro is defined for £ G ^°°(]R 3 , M 3 ) 
an ^ 

— [r ,£]u = (G(ro)f ■ R(r )^) - D(r )£ • R(r )u 

i=i 

where u = (ui, «2, U3) and £ = (£ , £ 2 , £ 3 ). 

PROOF. Let u G H s+1 (r,M 3 ). By definition of the surface rotational we have 

J R(r )u = -Trace(R(r )u). 
We deduce the 'rf 00 differentiability of i? and the first derivative in ro in the direction £ is 

— [r ,£]u = - Trace I — [r ,£]u 

= - Trace ( T [G(r )£] [R(r )u] ) - £>(r )£ • Trace (-R(r )u]) 



]T (G(r )& • R(ro)ui) - D(r )£ • i?(r )u. 



Here again we can obtain higher order derivatives of these operators iteratively. 

Remark 5.9 One can see that we do not need more than the first derivative of the defor- 
mations £. Thus these results hold true for boundaries and deformations of class c /f k+l , 
k E N* with differential operators considered in ^f c ('^ 7fc+1 (r), ( W k (T)). 
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(5.9) 



When ( = nwe obtain the commutators 

A(V r n)-V r (A.)=-^ r V r « 
d ( d \ 

— — (curlr u) — curlr — u = IZr curlr u — Hr curlr u 
On \on J 

^(divru) - divr U^uj = - Trace (ft r [V r u]) 
— (curlr u) — curlr ( ^-u j = — Trace (TZr [curlr u]) — Ht curlr u 

where Rr = [Vrn] and 7ir = divr n. 
From the precedent results we have: 

Lemma 5.10 The application 

L: -> ^ C {H S+2 (T),H S (Y)) 

r h-> TrAr.V 1 

is c to°° -Gateaux differentiable. 

PROOF. It suffice to write: 

T r Ar r Tr l = (r r divr r ^"^(rj-Vr^r" 1 ) = -(r r curlr, r" 1 )^ curl r? . r^ 1 ). 

The operators r r Ar r r~ 1 is composed of operators infinitely Gateaux differentiable. ■ 

View the integral representations of the operators PrC^P^ 1 and PrM^P" 1 we have 
to study the Gateaux differentiability of the applications r i— >■ r r Ap^ 1 curlr r t~ 1 and 
r i y Tr-Ap^ 1 curlr r r^ 1 . We have seen that for r G B^° the operator curlr r is linear 
and continuous from H s+1 (r r ) to H£(T r ), that the operator divr r is linear and contin- 
uous from TH s+1 (r r ) in H°{T r ) and that A^ 1 is defined from H*{T r ) in H s+2 (T r )/m. 
To use the chain rules, it is important to construct derivatives in r = between the 
spaces H s+1 (r) and H%(T) for the scalar curl operator, between the spaces TH s+1 (r) and 
H£(r) for the divergence operator and between the spaces H*(F) and H S+2 (T)/M. for the 
Laplace-Beltrami operator. As an alternative we use the : 

Proposition 5.11 Let u be a scalar function defined on F r . Then u r £ H^(T r ) if and 
only if J r T r u r = J r u r o (I + r) € H%(T). 

As a consequence the applications r i— >■ J r r r curlr r r^ 1 and r i-> J r r r divr r 7r _1 (r) are 
well-defined from H s+1 (r) and TH s+1 (r) respectively to tf* s (r). 

Lemma 5.12 The applications 

Bf -> ^ c (H s+1 (r),Hj(r)) ^ ^ c (TH s+1 (r),H:(r)) 

r J r r r curlr r r r 1 r (->■ J r r r divr r vr 1 (r) 
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are ff 00 - Gateaux differentiate end their first derivatives at r = defined for £ G < ^ 6 °°^ , 

<9J r r r curlr r r" 1 ^ 

— [0, f ]u = - 2^ Vr& • curlr ii*. 

i=l 



and 

dJ r r r divp r 7r _1 (r) 



[0,£]u = -Trace([V r £]V r u) + div r £div r u + ([V r u]n • [V r £]n) 



dr 

+ (u- [V r e]n)^r- 



PROOF. Let u G TH s+1 (r) We have 

<9r r 7r _1 (r) 



[0,£]u = (u-[V r £]n)n. 



<9r 

Next we use the lemma PTTI I5.5l and l5.8l For u G H S (T), it is clear that Yli=i Vr&curlr Hi 
is of vanishing mean value since the space Vr-ff s (T) is orthogonal to curlr H S (T) for the 
L 2 scalar product. An other argument whithout using the explicit form of the derivatives 
is : for all u G H S+1 (T), we derive the application 

r i — y J J r T r curlr r r 7 T 1 u da = 0. 

and for u G TH s+1 (r) we derive the application 

r i y J J r r r divrv 7r _1 (r)u da = 0. 

■ 

Let us note that u r G H s (T r )/R if and only if T r u r G H S (T)/R. 
Lemma 5.13 The application 

^ c (^ s+2 (r),^(r)/]R) 

is ff 00 -Gateaux differentiate. 

PROOF. We have seen in section [2j that the Laplace-Beltrami operator is invertible from 
H s+2 (r r )/R to H*(F r ). As a consequence J r T r Ar,,^ 1 is invertible from H S+2 (T)/R to 
i?*(r). By lemma 1331 we deduce that r i— > 7vAp r" 1 J^ 1 is ^"^-Gateaux differentiable 
and that we have 



Now we have all the tools to establish the differentiability properties of the electro- 
magnetic boundary integral operators and then of the solution to the dielectric scattering 
problem. 
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5.3 Shape derivatives of the solution to the dielectric problem 

For more simplicity in the writing we use the following notations : 

*E K (r) = * t eP;\ ^M.(r) = ^P- 1 , C K (r) = P r C r K P;\ et M K (r) = P r M r K P;\ 

Theorem 5.14 The applications 

B™ -> ^ c (TH-i(div r ,r),H(curl,K p )) 

r ^ ^>E K (r) 

are infinitely Gateaux differentiate. Moreover, their first derivative at r = can be ex- 
tended in linear an bounded operators from TH 2 (divp, T) to H(curl, $7) and H; oc (curl, f2 c ) 

1 d^E . m 

and given] € TH 2 (divp, T) the potentials — — — ^ [0, £lj et — - — - [0, £]j satisfy the Maxwell's 

or or 

equations 

curl curl u — k 2 u = 

in Q and fi c and the Silver-Muller condition. 

PROOF. Let j e TH~2(divr,r) and Vrp + curlrg its Helmholtz decomposition. Recall 
that ^E K (r)} and ^M K { r )\ can be written: 

=^r- 1 (r r P r - 1 j) - «- 1 V*;r- 1 (r r Ar r r- 1 p) ) 
^(r)j=curlCr- 1 (r r P r - 1 j). 

By composition of differentiable applications, we deduce that r 1— >■ *S>E K ( r ) an d r *— > ^> m k 
are infinitely Gateaux differentiable far from the boundary and that their first derivatives 
are continuous from TH2(divr,T) to L 2 (f]) U Lf oc (£l c ). Recall that we have, 

curl (r)j = K^ M A r )i and curl^M K (r)j = rc*£ K (r)j. 

Far from the boundary we can invert the differentiation with respect to x and the deriva- 
tion with respect to r and it gives: 

cur.^[0, 5 |j = ^ M et curl^O^j = ,*» M . 
d^>E 9^>m 

It follows that K 0, £ j and — — 11 0, £ j are in H(curl, fl) U H ioc (curl, Q, c ) and that 
Or Or 

they satisfy the Maxwell equations and the Silver-Muller condition. ■ 
We recall that the operator C K (r) admit the following representation : 

C K (r) = P r C£P r x j = V r P(r) + curl r Q(r), (5.10) 
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where 

P(r) = k (J^Ar.r- 1 )-^^ curler- 1 ) (tvVJt" 1 ) [(Ty-V^T" 1 ?) + (r r curlrv r^q)] 
and Q(r) = 

-K (J r T r Ar r T- 1 )- 1 (J r T r divr r 7r- 1 (r))7r(r) (tv^t" 1 ) [(tv-V^t" 1 ^ + (r r curlr r r^q)] 

Remark 5.15 Let] G TH s (divr,r) and Vr p + curlr g its helmholtz decomposition. 
We want to derive: 

PrClP; 1 ] = P r CX(Vr r T-V + curlr r r^q) 
= P r (Vr r P r + curlr r Q r ) 
= V r P(r) +curl r Q(r). 

We have: 

The derivative with respect to r a r de P r C£P r _1 j is given by the derivatives of the 
functions P(r) = T r (P r ) and of Q(r) = T r (Q r ). 
We also have 

By composition of infinite differentiable applications we obtain the 
Theorem 5.16 The application: 



B°° -> ££ r TH-2 



5(div r ,r),TH-5(div r ,r)) 



r i y PrClP' 1 

is infinitely Gateaux differentiable. 

Recall that the operator P r M£P r _1 admit the following representation : 

PrM^Pr" 1 ] = VrP'(r) + curlr Q'(r), 

where 

P'(r) = {JrTrA^r- 1 )- 1 («V r r r n r ■ (t^t" 1 ) [(rrVr^p) + (r r curlr r r^q)] 

+ {j rTr A Vr T- l Y l {J r T r DlT- l ){T r A Vr T- l P ) 

and Q'(r) = 

(JrTrAr^y 1 (J r T r curlr r V 1 )^^ - D^r" 1 ) [(r.Vr.r-^) + (r r curlr r r" 1 ?)] 
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with 



TrBlP' 1 ] = r r | jT VG(«, | • -y r |) (n r ( • ) • (Vr^-^Cl/r)) <^(l/r) 

+ y VG(k, I • -y r \) (n r (-) • (curlr r r T 7 1 g)(y r )) dcr(y r )j 



Theorem 5.17 TTie application: 

Bf 3 -> (TH^(div r ,r),TI-H(divr,r 
r h->. P r M K P- x 

is infinitely Gateaux differentiable and the Gateaux derivatives have the same regularity 
than M K so that it is compact. 

PROOF. By composition of infinite differentiable applications it remains to prove the 
infinite Gateaux differentiability of the application 

B 5 -> i? c (TH-5(div r ,r),H5(r) 

r ^ TrB^P- 1 . 

The function (x, y — x)\-t VG(k, \x — y\) is pseudo- homogeneous of class 0. We then have 

3 

to prove that for any fixed (x,y) £ (T x T)* and any function p € H2(T) the Gateaux 
derivatives of 

r i-^ (r r n r )(x) • (r r Vr r r r _1 p) (y) 
behave as \x — y\ 2 when x — y tends to zero. To do so, either we write 

(r r n r )(x) • {T r ^T r T~ l p) (y) = ((r r n r )(x) - (r r n r )(y)) • (r r Vr r T f T 1 p) (y) 

or we use lemmas 14.21 and 12.11 ■ 



Theorem 5.18 Assume that : 

1) E inc G H[ oc (curl,M 3 ) and 

2) the applications 

Bf -> TH-2(div r ,r) 



r P, (n r x (curl E inc ) |r J 



are Gateaux differentiable at r = 0. T/ien i/ie application mapping r onto the solution 
S(t) = E(O r ) £ H(curl,S7) U H/ oc (curl, f2 c ) to i/te scattering problem the obstacle O r is 
Gateaux differentiable at r = 0. 
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PROOF. By composition of differentiable applications. We write for the exterior field 



dr 



[o,£] 



+ 



+ 



+ 



+ 



- g r [o, £1 - ^ gr [o, £}c - vn^M Ke — [o, i] ) j 



<9r 

dM Ki inc dC Ki j 



dr 



-c« ^— M ■ 



The condition 1) guarantees that the solution j G TH2 (divr,T) so that the first term in 
the right handside are in H; oc (curl, fi c ) and the second condition guarantees that the last 
two term is in H; oc (curl, Sl c ). Of the same for the interior field we write: 



dS l 
dr 



[<U] 



1 d^p d^Af 
p-gf-M-f^ (E s + E mc ) - -£*■[(>, (E s + E-) 

1 T dp rlk (^ S (r) + E 



<9r 



9r 

0P r7 £ (^(r) + E-) 

■ [o, ej - *M Kj ^ [o, £j 



The condition 1) guarantees that 7^ (E s + E inc ) and 7^ (E s + E inc ) are in TH^ (div r , T) 
so that the first two terms are in H (cur 1, fi) and the second condition guarantees that 
the last two term is in H(curl, O). ■ 

Theorem 5.19 The application mapping r to the far field pattern E°°(fl r ) G ~Y c to°°{S 2 ) 
of the solution to the scattering problem the obstacle Q r is ff 00 - Gateaux differentiable. 

5.4 Characterisation of the first derivative 

The following theorem give a caracterisation of the first Gateaux derivative of r 1— > S{r) 
in r = 0. 



Theorem 5.20 Under the hypothesis of theorem \5.18\ the first derivative at r = in the 
direction £ solve the following scattering problem : 



curlcurl^-[0,e]-S 2 ^ r [0,e]=0 
curl curl — [0, £] - k\— [0, £] = 



(5.11) 
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with the boundary conditions : 



n x — [0,C] - n x — — [0,£] = g D 



dr 



dr 



(5.12) 



where 



Hi 1 n x curl— [0,f] - ^ e lf1 x curl [°> = 9N, 



+ curl r (£ • n) n • (E* - E s - E inc ) , 



and 



gN = - (f . n ) n x — (^r 1 curl E* - n' 1 curl E s - n' 1 curl E inc ) 
+ curl r (£ n) n (^r 1 curl E* - ^ curl E s - /xT/ 1 curl E inc ) . 

and where — — [0, £1 satisfies the Silver-Miiller condition, 
or 

PROOF. We have shown in the previous paragraph that the potential operators and 
their Gateaux derivatives satisfy the Maxwell' equations and the Silver-Miiller condition. 
It remains to compute the boundary conditions. We could use the integral representation 
as Potthast did but it would need to write too long formula. For x G V we derive in r = 
the expression: 

n r (x + r{x)) x {S\r){x + r(x)) - £ s {r){x + r{x)) - E inc {x + r(x))) = 0. (5.13) 
It gives in the direction £: 

= ^rMto x ( E ^) - ES (^) - Emc W) 
ir[°^]W- V [(U](x) J 

+n x (e(x) • V(E i - E s - E inc )) . 

Qt n 

We recall that ^ r [0,£](x) = — [Vr£] n and we use 

Vn = V r «+ ( ^ j n. 

We obtain : 

n(x) x (^rPUK*) - -£rPU]0«0) = [Vr£] n x (E*(x) - E'(s) - £**(*)) 

-n x • V r (E*(a;) - E s (x) - E" tc (x))) 
-(£ ■ n)n x ^- (P(x) - E s (x) - E iriC (x)) . 
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Since the tangential component of E l — E s — E mc vanish we have: 
(£0) • V r (E*(x) - E s (x) - E znc (x))) = ([ T V r n] £) (n • - E s (x) - E inc {x))) 

and 

[V r £] n x (E*(x) - E s {x) - E inc {x)) = ([V r £] n) x n (E*(a;) - E s {x) - E inc {x)) ■ n. 

For regular surface we have Vpn = T Vrn and 

([V r f] n) x n - n x ( [ T V r n] £) = curl r (£ ■ n) . 

We deduce the first boundary conditions. The second boundary condition corresponds to 
the same computation with the magnetic fields. ■ 

Using the commutators (|5.9p one can verify that the solution of this problem is in 
Hj oc (curl, O) since the trace 

d 

u € H^curl.fi) H- n x — u € ThT5(div r ,r) 

is linear and continuous. 

Conclusion 

In this paper we have presented a complete shape differentiability analysis of the solu- 
tion to the dielectric scattering problem using the boundary integral equation approach. 
These results can be extended to many others electromagnetic boundary value problems. 
Thanks to the numerous computations of Gateaux derivatives we obtain two alternatives 
to compute the first shape derivative of the solution : either we derive the integral repre- 
sentation or we solve the new boundary value problem associated to the shape derivatives 
with boundary integral equation method. Whereas this last alternative needs boundaries 
of class ff 2 at least since it appears any derivatives of the normal vector, many results in 
this paper are still available for Lipschitz domains as for example the computations of the 
Gateaux derivatives of all the surface differential operators of order 1 with deformations 
of class c ta l only and other functionals viewed in section 2J One can find in the litterature, 
the theory of pseudo-differential operators on Lipschitz domain |34] . it remains to find 
the optimal regularity of the deformations in order that this integral operators are still 
Gateaux differentiable. According to the Helmholtz decomposition we have on Lipschitz 
domain: 

TH~5(div r ,r) = V r (ft(O)0 curlr (h*(T)/R \ . 

where 

H(T) = {ue JZ^rJVR; A r € H^(T)}. 

If we want to extend the result to Lipschitz domain, we have to construct another invert- 
ible operator between H(T r ) and ~H(T). 
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